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The B-meson distribution amplitude (DA) is defined as the matrix element of a quark-antiquark 
bilocal light-cone operator in the heavy-quark effective theory, corresponding to a long-distance 
component in the factorization formula for exclusive B-meson decays. The evolution equation for the 
B-meson DA is governed by the cusp anomalous dimension as well as the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi-type anomalous dimension, and these anomalous dimensions give the "quasilocal" 
kernel in the coordinate-space representation. We show that this evolution equation can be solved 
analytically in the coordinate-space, accomplishing the relevant Sudakov resummation at the next- 
to-leading logarithmic accuracy. The quasilocal nature leads to a quite simple form of our solution 
which determines the B-meson DA with a quark-antiquark light-cone separation t in terms of the 
Q\ , DA at a lower renormalization scale with smaller interquark separations zt (z < 1). This formula 

allows us to present rigorous calculation of the B-meson DA at the factorization scale ~ \J 771(,Aqcd 
for t less than ~ 1 GeV" 1 , using the recently obtained operator product expansion of the DA as the 
input at fi ~ 1 GeV. We also derive the master formula, which reexpresses the integrals of the DA 
at /i ~ 'mbAQCD for the factorization formula by the compact integrals of the DA at fi ~ 1 GeV. 

D , PACS numbers: 12.38. Cy, 12.39.Hg, 12. 39. St. 13.25.Hw 
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The B-meson light-cone distribution amplitude (LCDA) is one of the important ingredients of the QCD factoriza- 
tion formula for exclusive B decays [Il-Q and has recently attracted much attention due to its central role for the 
analysis of the experimental data, e.g., hadronic and radiative B-decay data The B-meson LCDA appears in 
the factorization formula for hard spectator interaction amplitudes, where a large momentum is transferred to the 
spectator light-quark via gluon exchange flHl"fl| , and represents the nonperturbative matrix element that describes the 
' leading amplitude to have the valence quark and antiquark with a light- like separation inside the B meson pd| . Grozin 
y—{ ', and Neubert [l2j studied constraints on the B-meson LCDA from the equations of motion, heavy-quark symmetry 
!* ■ and the renormalization group (RG), and they gave the first quantitative estimate of the LCDA using QCD sum rules 
. £h ! with the leading perturbative and nonperturbative effects taken into account. The light-cone QCD sum rules for the 
B-decay form factors were also used [H|,[l4| to estimate the first inverse moment [HI of the LCDA, which participates 
in the corresponding factorization formula. The Grozin-Neubert's QCD sum rule calculation was extended by Braun, 
Ivanov and Korchemsky [l6| including the perturbative and nonperturbative corrections, and the importance of the 
NLO perturbative corrections was emphasized. Indeed, the true non-analytic behavior of the B-meson LCDA associ- 
ated with the "cusp singularities" [I3l is only revealed at this level including the radiative corrections [H], and it is 
this behavior that renders the (non-negative) moments of the LCDA divergent, even after renormalization (l2j (see 
fl9j for a similar behavior in three-quark LCDAs for the Ah baryon). Introducing the regularization for the moments 
with an additional momentum cutoff, Lee and Neubert [2(3 ] evaluated the first two moments for a large value of the 
cutoff in terms of the operator product expansion (OPE) with the NLO perturbative corrections, as well as the power 
corrections which are generated by the local operators of dimension 4, and they used the results as constrains on 
model building of the B-meson LCDA. 

Another feature peculiar to the B-meson LCDA is that it involves a complicated mixture of the multiparticlc Fock 
states of higher- twist nature through nonperturbative quark-gluon interactions, as demonstrated using the equations of 
motion and heavy-quark symmetry [l2l. l2ll. I22I] . A first systematic treatment of the mixing of the multiparticle states, 
disentangling the singularities from the radiative corrections, has recently been accomplished by the present authors, 
and the B-meson LCDA is obtained in a form of the OPE as the short-distance expansion for the quark-antiquark 
light-cone separation, with the subleading and subsubleading power corrections, generated by the local operators of 
dimension d = 4 and 5, respectively, and the NLO corrections for the corresponding Wilson coefficients [23[. This OPE 
enables us to evaluate the B-meson LCDA for interquark distances t with t < where \x is the renormalization 
scale of the LCDA, in a rigorous way in terms of three nonperturbative parameters in the heavy-quark effective theory 
(HQET), one of which is the usual mass difference between the B-meson and 6-quark, A = mg — mp, associated with 
matrix elements of dimension-4 operators, and the other two are the novel HQET parameters [HI, [2l|, [22[ associated 
with matrix elements of the quark-antiquark-gluon three-body operators of dimension 5. Note that the range of t where 
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the OPE is directly applicable becomes wider for the smaller value of the scale fi, as t < 1/fi: choosing fj, = 1 GcV, 
corresponding to typical hadronic scale, the model- independent result for interquark distances t < 1 GeV -1 has been 
obtained from the OPE and this result has also been used to constrain the behavior of the LCDA at large distances 
t > 1 GeV -1 [23[ . Those results of the S-meson LCDA at /i = 1 GcV have to be evolved to the factorization scale of 
order /ih c ~ \/mb^QCD that corresponds to the characteristic "hard-collinear" scale for hard spectator scattering in 
exclusive B decays [l|43( , when we substitute the LCDA into the relevant factorization formula. 

For this purpose, in principle, we can utilize the analytic solution for the evolution equation of the B-meson LCDA 
obtained in [13, 120]. However, the corresponding solution is directly applicable when the LCDA is given in the 
momentum representation, which we find inconvenient in our case: the Fourier transformation of the above OPE- 
based results to the momentum space mixes up the model- independent behavior for t < 1 GcV" 1 with the behavior 
for t > 1 GcV" 1 which relies on a certain model for the large t behavior. On the other hand, it has been noted 
that the relevant evolution kernel embodies the particularly simple geometrical structure in the coordinate-space 
representation These facts motivate us to treat the evolution of the i?-meson LCDA in an unconventional 

way, working in the coordinate-space representation. We are able to find the analytic solution for the corresponding 
evolution equation, and demonstrate that the solution determines the B-meson LCDA in terms of the LCDA at a 
lower scale /i with smaller interquark separations and thus preserves the boundary at t ~ 1 GeV -1 between the model- 
independent and -dependent behaviors of our LCDA, even after evolving from fi = 1 GcV to /ihc- Wc emphasize that 
such simple RG structure of the £?-meson LCDA can be manifested only in the coordinate space. Furthermore, as 
we shall demonstrate, it is this simple structure that enables us to derive the master formula, by which the relevant 
integrals of the LCDA at the scale //he, arising in the factorization formula for the exclusive B-meson decays, can be 
reexpressed by the compact integrals of the LCDA at the scale n = 1 GeV. Therefore, we believe that the coordinate- 
space approach for the RG evolution of the i?-meson LCDA deserves detailed discussions in the present paper. We also 
show that our solution can be organized so as to include the Sudakov rcsummation to the next-to-lcading logarithmic 
(NLL) accuracy, taking into account the effects of the anomalous dimension at the two-loop level which is associated 
with the cusp singularity. We present the first rigorous result of the B-meson LCDA at the relevant factorization 
scale fih c for t < 1 GcV" 1 . Combining with the results for the long-distance behavior, we discuss an estimate for the 
inverse moments of the LCDA at fi = /ihc- 

The paper is organized as follows. Sec. |TT] is mainly introductory; we give the operator definition of the i?-meson 
LCDA, explain the result for its renormalization in the coordinate space, and derive the corresponding RG evolution 
equation. We demonstrate in Sec. IHII that. as the solution of this equation, we can obtain the new coordinate-space 
representation for the evolution of the i?-meson LCDA, which manifests the simple RG structure, and also organize 
the result so as to include the Sudakov rcsummation at the NLL-lcvcl. In Sec. IIV1 we use our coordinate-space 
representation of the evolution to derive a compact and closed formula for the inverse moments of the LCDA in terms 
of the certain integrals of the LCDA at a lower scale /1. Application of our results to calculate the evolution of the 
OPE-based form of the B-meson LCDA is presented in Sec. [V] and we discuss an estimate for the inverse moments 
of the LCDA. Sec. IVII is reserved for conclusions. 



II. DEFINITION AND RENORMALIZATION IN THE COORDINATE SPACE 



The leading quark-antiquark component of the £?-meson LCDA is defined as the vacuum-to-meson matrix element 
in the HQET [lj: 



1 



\q(tn)[tn,0]i~f 5 h v (0)\B(v)) 



(1) 



where q(tn) is the light- antiquark field, h v (0) is the effective heavy-quark field, and these fields form a gauge-invariant 
bilocal operator linked by a light-like Wilson line, 



[tn, 0}=V exp 



ig dX n ■ A(Xn) 



(2) 



with n M as the light-like vector, n 2 = and n ■ v = 1, and w M denoting the 4- velocity of the B meson. The bilocal 
operator is renormalized at the scale /1 and, here and below, /1 refers to the MS renormalization scale. In the 
definition (J]), 



F(fi) = -i(0\qi~/ 5 h v \B(v)) 



(3) 



denotes the B-meson decay constant in the HQET [24[ and 0+(u;,/i) in the RHS is the LCDA in the momentum 
representation where ujv + denotes the light-cone "+" -component of the momentum carried by the light antiquark. 
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FIG. 1: The Feynman diagrams relevant for the one- loop renormalization of the nonlocal light-cone operator in {T}. The 
dashed line represents the Wilson line in between the quark fields, and the double line represents the effective heavy-quark 
field. 



The renormalization of the bilocal operator of ([1) was studied in fl2l , [l8j , calculating the UV divergence in the 
one-loop diagrams of Fig. 1 in the momentu m sp ace (see also [1^, |26|). The calculation of those diagrams has also 
been carried out in the coordinate space [Till [23| , and the result yields the renormalization of the bilocal operator 
®(t) = q(tn)[tn, 0] ^75/1^(0) in the coordinate-space representation as (unless otherwise indicated, a s = a s (fi)) 

e--(t) = e~(t,rt + {(-i - § + i) e™(t,„) 

4r <is (T^) + e "" (rf -" ) }- (4) 

connecting the bare and renormalized operators by the "z-dependent" renormalization constant in D = 4 — 2e 
dimensions and Feynman gauge, where Cf = (iV| — l)/(2N c ), and 

L = ]n[i(t-iO)fie 1B ] , (5) 

with the Euler constant and the "— «0" prescription coming from the position of the pole in the relevant propagators 
in the coordinate space. The plus-distribution is defined, as usual, as 

'„(VM /M.f* M , (6) 



for a smooth test function f(z). In the one-loop contributions in (j4|, the first two terms, the double-pole term and 
the single-pole term involving L, manifest the cusp singularity [l?} in Fig. 1 (a), i.e., the singularity in the radiative 
correction around the cusp (at the origin) in the Wilson line, 

[in, 0] [0,-ooi;] , (7) 

which is contained in ([1]), using the relation h v (0) = [0, — oov]h v (— oov). The last term in ([J) comes from Fig. 1 (b), 
accompanying the plus-distribution characteristic of the loop integral associated with the massless degrees of freedom 
only, while the remaining one-loop term comes from the contribution of the renormalization constants of the two 
quark fields, q and h v . We note that Fig. 1 (c) is UV-finite in the Feynman gauge [27J and does not contribute to ([4]). 

The RG invariance to the one-loop accuracy based on (j4]) implies that the £?-meson LCDA ([I]) obeys the evolution 
equation in the coordinate space, 

d * 



^—<p+(t,iJ,) = -[r cnsp (a s )L + ^F(a s )](j} + (t,fj,)+l dzK(z,a s )4> + {zt,fi) , (8) 



with the one-loop RG functions 



r (a) = r (1) — r (1) - 4cv (9) 

1 cusp V"s; — L cusp a J L cusp — j \VJ 

lF(a s ) = ^^, 7 ^=2C F , (10) 



4~ 



and 



(11) 
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Here, r cusp (a s ) corresponds to the anomalous dimension of the Wilson line with a cusp, Q, and coincides with the 
LO term of the universal cusp anomalous dimension of Wilson loops with light-like segments [17| ; we obtain the finite 
result ©, because the contribution from the double-pole term of Q is canceled by the contribution generated by 
taking the derivative of L in the next term with respect to \i. 7f(a s ) of (|10[) represents the anomalous dimension 
from the above-mentioned contribution of the renormalization constants of the two quark fields, combined with the 
so-called hybrid anomalous dimension of heavy-light currents in the HQET [24| , which governs the scale dependence 
of the decay constant ((3]) as, at one-loop accuracy, 

^t/^ = 3Cf S f(m) ■ (12) 

K(z, a s ) of (jTTJ) comes from the last term of (@| and represents the z-dependent anomalous dimension associated with 
the massless degrees of freedom only. We note the remarkable property in © that the evolution kernel in the RHS, 
composed of (f9"))- (|ll[) . is quasilocal, such that the evolution mixes the LCDA with itself and with the LCDA associated 
with smaller light-cone separation zt (z < 1). This is due to the similar structure appearing in the renormalization in 
the coordinate space, ©, and reflects [16| the fact that the cusp renormalization induced by Fig. 1 (a) is multiplicative 
in the coordinate space [l7| while Fig. 1 (b) gives the contribution identical to the similar correction to the light- 
quark-antiquark bilocal operators, which embodies simple geometrical structure in the coordinate space so as to mix 
operators associated with smaller "size" [28| . 

It is straightforward to perform the Fourier transformation of (JSJ) to the momentum space and derive the evolution 
equation for 4>+(uj, M) of |T]), using 



1 

2^ 



dte wt L<j) + (t,n) 

-OO 

Z7T Jq \Ld' — UJ — l\J fj, UJ' — UJ + %{] /X / 

-</> + (w,/i)ln / duo' — [<j) + (w',ij) -4> + (ui,fj,)} , (13) 



/i J LO — LO' 

and 



27T ~' J Q ~~ V 1 — Z 



±- I dte iut I dz[-— I 4>+(zt,n) 



<M^,m)+ I du'- e{u) ' u) [M"',»)-M",»)] , (14) 



UJ' UJ' 



and, indeed, the result reproduces the evolution equation obtained through the renormalization of the bilocal oper- 
ator of |T]) in the momentum space HH, H|| . We note that the momentum representation of the kernel in (fT4")) 
coincides with (a part of) the Brodsky-Lepage kernel for the pion LCDA [29[ and physically represents a Dokshitzer- 
Gribov-Lipatov-Altarclli-Parisi (DGLAP) splitting function that vanishes for lo/cu' > 1. On the other hand, a local 
contribution in ([5]). associated with the cusp anomalous dimension r cusp (o; s ), yields the new evolution kernel for 
lo/uj' > 1 in the RHS of (|13[) . so that the evolution in the momentum space mixes the LCDA 4>+(oj, (J.) with <fi + (oj', fx) 
over the entire region, < lo/uj' < oo [H, Hll- We also note that we cannot derive the moment-space represen- 
tation of the evolution equation (JSJ in a usual way as in the case of the LCD As for the light mesons [3(3|, because 
the presence of the logarithm ([5]) prevents us from performing the Taylor expansion of flSJ about t — 0; indeed, (|3]) 
shows that the renormalized LCDA is non-analytic at t = (see also the discussion in [l6|,[23(). Thus, the evolution 
equation for the i?-meson LCDA manifests simple geometrical structure as the quasilocality of the kernel only in the 
coordinate-space representation l[5]). 

One may anticipate that the evolution equation ((SJ) would hold to all orders in perturbation theory by taking into 
account the higher-loop terms in the RG functions ([^- (fTTj) . This is correct, at least, for a particular class of higher- 
loop corrections associated with the universal cusp anomalous dimension r cusp (a s ) of Wilson loops. For example, 
when we take into account the dia gram s that correspond to the two- loop corrections to the relevant Wilson line (J7]), 



r c usp(o; s ) of © gets modified into [17 1 



Tcusp V 



( a )- r « ^i + r^ (^Y (15) 
with (a3 



r< 2 ) = 4cv 

cusp w r 



67 tt 2 \ _ 10 Ar 

C G N f 

9 3 J G 9 1 



(16) 
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where Cq = N c and Nf denotes the number of active flavors. Actually, it is not known at present whether the 
effects of all the other two-loop corrections to the bilocal operator in (TTJ) can be absorbed into the remaining two RG 
functions in ©, 7f(q ; s ) and K(z,a s ), as their two-loop terms. Still, the evolution equation ((HJ, with ([TU)l . (fTT|) and 
the two- loop cusp anomalous dimension (|15[) taken into account, is useful for resumming the Sudakov logarithms to 
a consistent accuracy, as we will demonstrate in the next section. 



III. ANALYTIC SOLUTION IN THE COORDINATE SPACE 



The LO solution for the evolution equation of the i?-meson LCDA was obtained in the momentum representation in 
[HI, and the result determines /x) of (fl} as the convolution of (f> + (oj' , fio) at a lower scale /io and the (complicated) 

evolution operator, over the entire range of ui' (see (|A7j) in Appendix A). Its Fourier transformation in principle gives 
the solution for our evolution equation, (|5])- (fTTj) . in the coordinate space, but, in practice, we find it more useful to 
solve © directly: mathematically, © is an integro-diffcrential equation of similar type as the corresponding equation 
in the momentum space and has simpler structure for the kernel of integral operator than the latter case, as noted in 
Sec. HU This would imply that the strategy devised in [l8| to solve the evolution equation for the latter case should 
also allow us to solve ([5]). possibly with simpler manipulations. Moreover, intermediate steps of those manipulations 
reveal peculiar structures behind a rather simple final form of our solution, (|35[) below. 

First of all, we demonstrate that the strategy of [l8| is applicable to ((8j) and allows us to construct its general solution 
which is exact even when the higher-loop terms in the RG functions r cusp (a s ), 7i?(a s ), K(z, a s ) are taken into account. 
For this purpose, we further put forward the above-mentioned similarity between ([5]) and the corresponding integro- 
diffcrential equation of [l8j in the momentum space, by performing the analytic continuation for ([5]) as t — > —it. Then 
the evolution equation ([8]) becomes the integro-differential equation for the B-meson LCDA at imaginary light-cone 
separation, 4>+{— it, fi), as 

d 



ix— 4> + {-it,h) = - [r cusp (a s )ln(r / ue 7B ) + -f F (a s )](j) + (-iT, /i) + / dzK(z,a B )<j>+(-iZT, fi) . (17) 
afJ> Jo 

We recall that the kernel K(z, a s ) corresponds to the coordinate-space representation of a DGLAP-type splitting 
function and thus can be diagonalized in the moment space. In the coordinate-space language [28| , the corresponding 
moment is given as 

K(j,a s )= dzz j K(z,a s )=IC^(j)-^+--- , (18) 
Jo 4?r 

where (fTTj) gives the coefficient for the order a s term as 

K<U(j) = AC F jf dzz j (jZT^) = - 4C fH'U + 2) + 7J5 - 1] , (19) 

with ip{z) = (d/dz) lnT(z) being the di-gamma function, and the ellipses in (|18[) stand for the (presently unknown) 
terms of order a 2 s and higher. As mentioned below (|14[) . however, the usual moment is not useful for treating (|17l) : 
the presence of logarithm ln(r/ie 7i5 ) in the RHS suggests that the values for the moment j will be modified under the 
variation of the scale pi. The authors in [l8[ demonstrated that taking into account the corresponding "evolution" of 
the moment j indeed enables them to construct the general solution of the corresponding integro-differential equation 
in the momentum space (see also [H|). Thus, we take the ansatz, 

-i />c+ioo 

M-ir,v) = —. dj (rfx e^y- (M <p(j,Li) , (20) 

27T« Jc-too 

with a real constant c, and we determine <p(j,/j,) and £(A*,Ato) such that (|20[) obeys (|17[) . This ansatz has the form 
similar to the inverse Mellin transformation to construct the solution for the DGLAP-type evolution equation in the 
coordinate-space language (see [28j|). except for the contribution of £(/i, /io), which describes the evolution of the 
power of t from a certain (low) scale po to the scale p. We assume £(Mo, /io) = 0, without loss of generality, and /iq 
multiplied by e lE is put in the integrand of (|2"0"]l for convenience. Substituting (|2T)|) into (fT7|) . we obtain 



a/j, 



r c usp(a s )hi(i7/e 7E ) - j F (a s ) +fC(j - C(i",A*o), a s ) 



A* 1 ln(r/i e 7E t 



V>U,H) ■ (21) 
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Because the RHS of this equation should be independent of r, /io) obeys 



rcusp(^s) 



(22) 



This shows that /xq) is independent of j and is integrated to give, introducing the /3 function, /3(a s ) = /ida s /d/i, 



».(c) 



cusp 



(a) 



da = S(a s (/x), a s (/Jo)) 



Now (|2T1) reduces to 



dfj, 



-r cusp (a s ) In 7F(a s ) + K(j- Mo), a s ) 

Mo 



¥>(j, m) 



and this simple differential equation is immediately solved to give 

<Pti> /") = exp [V(/i, Mo) + W(/i, Mo, i)] <^(j, Mo) 

where 



<x.(mo) /%) 



r 



cusp 



(a) 



da 



a.(«0 /A" ) 



;(PC>) 



/3(a) 



^(j-S(a,a s (/i )), a) 



(23) 

(24) 

(25) 

(26) 
(27) 



and <p{j,Ho) should be expressed by the Mellin transform of the initial condition, </>+(— it, (j,q), from ([2TH) (see ([2^)1 
below). Substituting these results into (f20|) . we obtain 



+ (-iT,/i) = e v ^^ o) (^oe 7B )" ? 



— 0+(-«r ,/i ) 



2?ri 



,W(/i,/i ,j) 



(28) 



Here and below, £ = £(/x,/zo), unless otherwise indicated. The formula (|28[1 in principle gives the solution for 
(|17p. which is exact even when the higher-order terms in r cusp (a s ), 7j?(a s ), iT(M,a s ) are taken into account. How- 
ever, this solution has been obtained by assuming tacitly that (p(j,fi) of (f20|). expressed by the Mellin transform of 
(r/j e 7E ) 4 4>+(-iT,n) as 



= / — (rxtoe 7 - 5 ) J (Tii e 7E )^ + (-ir,/i) , 



(29) 



is a regular function in a certain "band" of region in the complex j plane, and that the constant c in (|28P is chosen 
such that the integration contour in this formula is contained within this band. Now we consider the condition for the 
convergence of the integral in (|29|). which in turn determines this band, as well as the range where (|28| is applicable: 
the short-distance behavior of 4> + {—iT, fi) as r — > in the integrand of (|2"9"]l can be determined [2(J[23| by perturbation 
theory, as a constant modulo lnr (sec ([65"T) below), so that ([29| is convergent for the integration region r ~ when 
£ > Re(j'). On the other hand, studies of the IR structure of the DA indicate </>+(— it, fi) ~ 1/r 2 or more strongly 
suppressed as r — > 00 [H, [2l|, so that the integral in (|29|) is convergent as r -» 00 when £ < Re(j') + 2. These 
considerations show that (|2TJ)) gives a regular function for the band with £ — 2 < Rc(j) < £ in the complex j plane, 
and the constant c in (l20l) . (|28l) should be chosen as 



£-2<c<£ . 



(30) 



Because £ of ([23| grows from 0, as \x increases from fio (see (0, (fl~5j) . and (fl6|)). only for the values of scales /x and Mo 
satisfying 

C = £(M,Mo)<2, (31) 

c can be chosen as a fixed constant and thus the solution ([28| describes the exact evolution of the B meson LCDA 
from fio to /i. (Note that the condition for the convergence of the convolution integrals of (|28|) and the corresponding 
hard part in the QCD factorization formula for exclusive B decays eventually requires (|60p below.) 
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To proceed further, we change the integration variable in (|27[) from a to H(o!, a s {fio))- Defining a x such that 
£ (a x , a s (//o)) = a;, we obtain 

HWo,j)= fdx ^f^ °! x) = [' dx ^^-^ +•■■ , (32) 

where the ellipses stand for the NLO or higher contributions that involve the two- or higher-loop anomalous dimensions. 
Using (0 and (19]), one finds 

e " e r(i + 2) ' (, " J 

up to the corrections of the two- loop level. In the complex j plane, (|33|) has poles at j = £ — 2 — n with n = 0, 1, ■ • • , 
which are all located in the left of the integration contour in ([28]) with (|30|) : by the theorem of residues, these poles 
give rise to nonzero contribution to the j integral for r > r' , while the j integral vanishes for t < r'. Evaluation of 
those pole contributions yields 



c+ioo 



d0_ / Ty r(j + 2-Q = _ (-1)" (T^-Z-n 

2m\T') T(j + 2) 1 ' ^ n\T(£_~n) W 



= 0(r-r') - ^Z-,^ ■ (34) 



(l-r'/r) 1 -^^) 

Substituting this result into ([28]) and changing the integration variable from t' to z = t'/t, we obtain 

p (i-tb)£ /■! / z \ 1_ £ 
+ (-iT,M) = e V( ^ o) (r/zoe^)^ ^ ^ J + (-i«,po) , (35) 

which is exact up to the NLO corrections mentioned in (|32|) . The contribution of the kernel (fTTj) in perturbation 
theory receives the RG improvement in (|35|) as (z/[l — z]) 1 ^^ with the modified power 1—0 where £ of ([23)) is 
induced by the cusp anomalous dimension. For the case with £ — > 0, we have the singular behavior as 1/(1 — z) 1 ^ = 
(1/£)£(1 — z ) + 1/(1 — z)+ + O(0, but this eventually gives the finite contribution to the RHS of (|35l) . combined 
with the £ — >• behavior of the gamma function, T(0 = l/£ — Je + 0(0- This also shows that the RHS of (|3"5)) 
reduces to <p + (—ir, fi ) when // — > i-e., when £ — >• and V(/x,/x ) — > (see (f23|) . (|26|) ). as it should be. In (f35|) . 
it is straightforward to perform the analytic continuation from the imaginary light-cone separation to the real one, 
as r — > it, and the resulting solution for the evolution equation (JSJ embodies a quite simple structure to determine 
the -B-meson LCDA with a quark-antiquark light-cone separation t in terms of the LCDA at a lower renormalization 
scale [1q with smaller interquark separations. The Fourier transformation of this result is calculated in Appendix fAI 
and the obtained momentum representation (|A7[) reproduces the result of [IH [2p| derived in the momentum space; 
in particular, the factor (r/ioe 7E )~^ in (|35[) , which is non-analytic at r — > 0, produces the radiative tail as ~ uj^^ 1 
for large in (|A7[) . which renders all non-negative moments of the LCDA, J °° duiuj n (f>+(uj, /x) with n = 0,1,2,..., 
divergent, irrespective of the initial behavior, <^+(w , //o) [HI- We a l so emphasize that our result (|35[) has a much 
simpler structure than (|A7[) ; i.e., the most compact expression possible for calculating the evolution of the B- meson 
LCDA under changes of the renormalization scale is provided by our coordinate-space result (I35|) . 
We note that the first two factors in (|31)|) . given by 

e V(p,m>) ( TMoe T B )-« = e V(M,Mo)-Sln(T Mo e^) ^ ^ 



are unaffected by the above manipulations (|33|) . (|34|) . which are valid up to the NLO corrections, and thus (|36|) gives 
the exact result even when the higher-order terms in r cusp (a s ) and 7^ (a s ) are taken into account for (|23p and (|26l) . 
Indeed, substituting those definitions of £ and V(^t,/xo), we may reexpress the exponent of the RHS in (|36p as 

V( /1 , M o)-^ln(w^) = - -^[r cusp (a s ( M '))ln(V e ^)+ 7 F(a s (M'))] , (37) 

which shows that ([55]) corresponds to the general solution of the evolution equation (fTT)) with the contribution of 
the kernel K(z,a s ) omitted |32|]. We now derive the explicit form of (|2"3"|) and (|2^|) arising in our solution (|2"8"|) . (|35[) : 
substituting (j9)), (jTUj) . (|15|) and the usual perturbative expansion for the /3 function, 

«*>-*£--*-i>(2r- 



^ „=0 

3 3 ' ' M 3 G 3 



A) = — C G - -N f , = -C% - — Cety - 20^ , • ■ ■ , (38) 
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a straightforward calculation gives 



£(m.mo) 




(39) 



and 



V(//,/x ) 




a s (//) 




(40) 



where the ellipses stand for the terms that are down by a s compared with the preceding terms and receive the 

(3) 

contributions due to higher loops, e.g., the three- loop cusp anomalous dimension Tcuspj the two- loop local anomalous 

(2) 

dimension j F , etc. If we substitute only the one-loop terms of these results, the first term of (f3"9")l and the first line 
of (|40|) . into (|35|) . we obtain the explicit analytic form of the solution for the evolution equation (|SJ), exact at the 
one- loop level with (|9|)- (fTTj) . 

The definition (ffij)) shows that V(fJ,, /io) involves the contribution associated with the first term in the RHS of 
(|24[) . i.e., the cusp anomalous dimension accompanying ln(/i//io) ~ l/ct s . As a result, in (|40p. the contributions 
associated with the cusp anomalous dimension are enhanced by the factor induced by this logarithm, compared to 
the contribution from the second term of (|2"rj)) with the local anomalous dimension (flU)) : the leading term is given by 

the one- loop cusp anomalous dimension Tculp, while the one- loop local anomalous dimension jp contributes to the 
next-to- leading term, i.e., at the same level as the two-loop cusp anomalous dimension. Here, the contribution due 
to 7L 1 ' corresponds to the one- loop level in the usual RG-improved perturbation theory, and thus the treatment at 
this level has to be complemented with the two-loop contributions associated with the cusp anomalous dimension, the 
second line of (|4"0)) . This pattern is characteristic of the Sudakov-type large logarithmic effects induced by the cusp 
anomalous dimension. This fact also requires us to reorganize our result (|35[) as well as (|40[) according to consistent 
order counting of those logarithmic contributions. This can be achieved by introducing 



a g (/z) 

X = Po—. In — 

4tt 



(41) 



and by following the standard procedure used in the soft gluon resummation formalism in QCD [33||: we organize (|4"0)) 
by a systematic large logarithmic expansion, where x ls formally considered of order unity and the small expansion 
parameter is a s (/j,), leading to 



In 



ln(l - x) - 



Substituting this expansion, (|40|) is recast into 



a s (n) Pi ln(l - x) 
4tt fa 1-x 



+ 0(a 2 s 



(42) 



V(//,Mo) 



47T 



(43) 



up to the corrections of 0(a s ), with 



h w (x) 



r 



(i) 



^[(X-I)ln(l- X )-X] 



r (i) 

1 cusp 



2ft 



ln 2 (l-x) + 



.(2) 
cusp 



cusp 



(ln(l -X)+X) 



+ I^ln(l- X ) 



(44) 
(45) 



In (|4"3"j) . the first and the second terms, (An/a s )h^ (x) and h^^x), collect the terms a™ ln n+1 (p 2 / fj.^) and 
a™ hi n (/j, 2 / /Iq) , respectively, with n = 1,2,..., corresponding to the LL and NLL contributions, while the 0(a s ) 
corrections omitted from (|4"3")l correspond to the NNLL or higher level. Thus, the first factor e v ^ AI,AI0 - ) in (|3"5|) is the 
exponentiation of the logarithmic terms a" ln m ( u? J n 2 ) with m < n + 1, playing analogous role as the Sudakov form 
factor in the soft gluon resummation in QCD [33|. It is straightforward to see that this factor e v ^^ ^ with only 
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the LL term, (A-k / a s )h^ Cl \x) = — (a s /327r)rcu sp 1ii 2 (/^ 2 //^q) +•■•, retained in the exponent (|43| corresponds to the 
double leading logarithmic approximation summing up the towers of logarithms a™ In "(/i 2 /^ 2 ,), and e v ^^ ' with the 
exponent (|43| at the NLL accuracy resums the first three towers of logarithms, a™ ln m (/i 2 //x§) with m = 2n, 2n — 1, 
and 2n — 2, exactly, to all orders in a s . 

The logarithmic expansion (|4"2l can be also applied to (j3"9")l . yielding 

^,/io) = -^ln(l- X ), (46) 

up to the corrections of 0(a s ). This result does not receive the logarithmic enhancement but obeys order counting 
similar as the contribution from the second term in (|26p due to the local anomalous dimension jp, as apparent 
comparing (j2"3")l with (f2l))) . Thus, the substitution of P5|) into ([55]) produces the NLL- level contributions while the 
omitted 0(a s ) contributions correspond to the NNLL or higher level, using the order counting similar as in (|43|) . 
Therefore, our solution (|35|) with (|43|) - (|46p substituted embodies the evolution of the B-meson LCDA, accomplishing 
the relevant Sudakov rcsummation, and is exact up to the corrections of the NNLL-level. We note that controlling 
the NNLL-level effects completely requires to take into account the three- loop cusp anomalous dimension Lcusp, as 
well as the local anomalous dimension and DGLAP-type splitting function at the two- loop level, "fp and K.( 2 \j), in 

dH) with (ED), dnj). 

Before ending this section, we mention that the factor ln(r/xoe 7E ) accompanying £ in the RHS of (|36[) could produce 
an additional logarithmic enhancement. Also, in the integrand of (f3"5"j) . the behavior as z — > and z — > 1 could receive 
another logarithmic enhancement. These facts suggest that the higher-order terms associated with the similar types 
of logarithms could be relevant if we intend to determine the precise shape of the LCDA at the "edge" . However, we 
do not go into the details of such higher-order effects here: systematic treatment of those higher-order effects would 
require an approach, which is beyond the scope of this work based on the evolution equation for the renormalization 
scale /1. Furthermore, it is the integrals of the LCDA over r, like (|4"8")l below, that is eventually relevant to exclusive 
B decays, and the above types of higher-order logarithmic effects at the edge region should play minor roles on the 
value of those convergent integrals, where the only relevant logarithm is \n(p? / ^) treated in this section. 

IV. MASTER FORMULA FOR THE INVERSE MOMENTS OF THE LCDA 

The B-meson LCDA (TTJ) with [i = /j,h c participates in the QCD factorization formula for exclusive B decays (Tl. [8l-FTol| 
through the inverse moments, 

^sHm) = / — 0+( w ,Ai) > a n (fi) = Xb(p-) [ —(f>+{tJ,fj,)]n n —. (47) 

Here, A^ 1 (/i) appears as the convolution with the hard part in the LO for the hard spectator amplitudes, while the 
calculation of the NLO effects for the hard spectator amplitudes requires also the logarithmic moments a n (/x) with 
n = 1, 2. We introduce the logarithmic moments in the coordinate space, 







R n {^)= I dr^+(-ir,/x)ln"(r/xe 7B ) , (48) 

which arc related to (l47l) as 



A^O*) = -Ro(m) , tn(p) = X B (p)Ri(p) , T 2 (fj,) = \b{ix)R 2 (ij,) . (49) 
These relations can be obtained, e.g., by considering the generating function for the inverse moments (|47p . 

^y<M^HEf#r, (50) 

and relating this to the generating function for i? n (/i) of (|4"8")l . i.e., / °° dr(r/ie 7B ) s 0+(— ir, fi), as 

— fH) (j)+ (u J ^) = ——- dTW>y<t> + (-iT t ri, (51) 

<jj \uJ r(i + s) j Q 
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where e S7E /r(l + s) = cxp[— ^^/ 2 (— s) fc £(fc)/fc], with ((k) being the Riemann zeta-function. Remarkably, the 
simple form of our solution (|35[) allows us to express the generating function in the RHS of (|51j) as 



x (— ) [°° dT(TiM>e^) s -t4> + (-iT,» ) 
VMo/ Jo 



(52) 



which plays role of the master formula to derive all the relevant moments at fi = /ihc in terms of the integrals of the 
LCDA at a lower scale [1q: Taylor expanding the both sides of this formula about s = 0, we immediately find 



A B 1 ( M )=i? (/i) = 



3 V(m,mo) + (1-7-b)C f°° 



r(i + o 



dT(T(j, e~> E ) e c/> + (-iT,^ ) 



and 



p V(/x,Mo) + (1-7e)« f°° 



r(i + o 

; V(m,Mo) + (1-7e)? 

r(i + 



dr [\n(T^ E ) +V(1 +0] (r/ioe 7E )-^+Hr,/i ) 



x / dri [ln(T/. e 2 ^)+^(l + e)] 2 -^(l+0 + y |(TMoe 7B )-^ + (-iT,Mo), 



(53) 



(54) 



(55) 



where ip'(z) = (d/dz)ip(z); it is also possible to derive the similar formulae for R n (fJ-) {n > 3). These formulae 
(|53| - ([55|l . combined with (j49]l. allow us to calculate (|47|) with the LCDA </>+(— ir, //o) as the input at the hadronic 
scale, and the results are exact up to the NNLL corrections when substituting (j43|) -(j46 ]) . Furthermore, it is worth 
presenting the corresponding results transformed into the momentum representation: substituting (|52l) into (|51[) . we 
obtain our master formula in the momentum representation as 



dui / /i 



(-) <M<^/f) = e 



v( K ,»o)+fi-2^ E )£ r ( 1 ~ ^) r (! + a - fl 

r(i + s)r(i- s + o 



\HoJ J U \r- 



(56) 



and, using (f50|) . we obtain 



a* 1 00 

A B (/i) 



3 V(m,Mo) + (1-2 7 e)€ 



_ e V(M,Mo)+(l-2 7 J5)€ 



In 



27E 



3 V(/i,/i )+(l-27E)€ 



In 



/f.e 



27E 



r(i-o 


r 00 




r(i + f) 




U! ' 


r(i-0 






r(i + o 






■ + 


+ 


-^(1 


r(i-o 


I* °° 




r(i + o 


7o 


w 
-, 2 



+ V(i + + 



<M^, Mo) 



+ V>'(i-0-^(i + 



— <M^,Mo) , 
Mo/ 



(57) 



(58) 



(59) 



and so on. We note that (|57|) and ([58)) reproduce the corresponding results that were found in [34| in the context of 
the RG evolution at the one-loop level in the momentum representation, while the closed form (|59j) for 02, as well 
as the above formulae ([53 |) - ([55"T) in the coordinate-space representation, is new. We also emphasize that our master 
formulae (|52"j) . ([55]) allow us to derive the closed form for the higher logarithmic moments straightforwardly. 

The above results (|52|) - ([59|) demonstrate that, as a result of the evolution, the relevant logarithmic-moment integrals 
have the common, additional "weight functions" determined by £, i.e., (T^oe 7E ) - ^ and (w/no)^ m the coordinate and 
momentum representations, respectively, compared with the corresponding formulae for [i = /iQ. In particular, the 
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condition for the convergence of the integrals in (1531 - ([55)) . taking into account the r — > (and r 
4> + (—iT, fi ) mentioned above ([3"0")) . indicates that only for the values of scales m and Mo satisfying 



£ = £(m,mo) < i 



oo) behavior of 



(60) 



our formulae ([53)) -([55)) (and (|57)) - (|59|) 

M = Mhc ( 



59])) are well-defined and applicable. This is indeed satisfied for the relevant scales, 
yj rrib Aqqd) and Mo ~ 1 GeV, as discussed in the applications in the next section (see Fig. [4] below), and, 



actually, even for quite large values of \x. 

The perturbative expansion of the above results (f51?]) - (f55|) (or J5_ 



55)1 ) in terms of o: s (m) (= «s) yields (see (|1T)) . 



^bV) = A^Mo) 
^i(m) = A^^Mo) 

+#i(Mo 
iM#) = A^^Mo) 

+i?i(Mo 

+ i?2(M0 



1 - 



a s CF 


In 2 






Mo 


- -Ri(Mo) 






Mo 




7T 




In 2 


JL 




ln^ 

Mo 




Mo 


2tt 


Mo 


2tt 


+ 6 


2tt 


In 2 


Mo 


2tt 


In — 

Mo. 


- #2(Mo) 


7T 



ln^- 

Mo 



Mo 



(61) 



(62) 



M a s C F 3 m , Oi s C F 2 m , 7ra< s Ci? M , 2a s Cj? A/oN 
m In 1 In 1 In 1 £(3) 



Mo 



2tt 



Mo 



2 71 



Mo 



■In" 



Mo 



27T Mo 



a s C F 


In 


_M_ 


+ " 


7T 


Mo 




a s C F 


In 


_M_" 




2?r 


A*o_ 





3 mo 
In — 



In — 

Mo 



Mo 

Q; s Ci? m 

«3(Mo) In — 

7T Mo 



(63) 



etc., to the corrections of order a 2 ; here we have substituted ©, (|TU)) for the relevant anomalous dimensions. These 
relations can also be obtained by substituting the expansion (|68)l . discussed below, into the RHS of (|48)l . Note that 
the correction terms in (|6T)l - (|63)l have the relative size < (o^Cf/tt) ln 2 (^/^ ), compared to the leading terms. In 
particular, our first relation ([6T)l has the double logarithmic term, A^ 1 (mo)[— {cx s Cf ln 2 (/////o)], which was absent 
in the corresponding relation discussed in (l6| ; this term comes from the expansion of the LL term of V(m, Mo) (see 
(|43)) . (04)) ) and it is straightforward to check that, taking into account this double logarithmic term, (|6T)) satisfies the 
correct evolution equation for A^ 1 (m), which is obtained by integrating (fT7| over r (see also ([68)) below). The formulae 
(|ST)) - (|63p allow us to relate the relevant logarithmic moments at m = Mhc (see P7 )) - P9^) ) to the similar moments at 
Mo ~ 1 GeV. As we demonstrate in the next section, the formulae (|6T)) - (|63| actually show good accuracy at the 
relevant scales and thus provide model-independent relations that are expected to be useful for the phenomcnological 
applications. The pattern illustrated in those formulae is general, such that the calculation oi R n (fi) at the fixed-order 
a s requires the knowledge of the n + 2 logarithmic moments at moj -Rfc(Mo) (k < n + 1). 



V. EVOLUTION OF THE OPE-BASED LCDA 



In this section we apply the evolution represented by our solution (|35p to a suitable input LCDA at the initial scale 
Mo, and study the behavior of the resultant LCDA <j)+(t, fx) at a higher scale fi to clarify the effects of the NLL-level 
evolution quantitatively. For the input LCDA, we note that the model-independent information is now available 
based on the OPE for the bilocal operator in ([T]) as the short-distance expansion for the quark-antiquark light-cone 
separation t [23[, i.e., 

q(tn)[tn,0W l5 h v {0) ~ ^ d(t, M)Oi(M) , (64) 

i 

as t —> 0, with the local composite operators C^(m) and the corresponding Wilson coefficients Ci(t,fj,), depending 
on the (MS) renormalization scale (i for the bilocal operator. Apparently, the lowest-dimensional operator that 
participates in the RHS is given by the dimcnsion-3 operator appearing in the definition ^ for the decay constant 
F(fj,), and, here, a complete set of local operators of dimension d = 4 and 5 is also taken into account. The dimension 
counting tells us that the Wilson coefficients associated with the dimension-d operators behave as Ci(t,fx) ~ t d ~ 3 
modulo logarithm. Those coefficient functions are calculated to the NLO (0(a s )) accuracy, and the corresponding 
NLO corrections turn out to induce the contributions associated with the logarithm L of ([5j) . Substituting the result 
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for ([M]) into ([T]), the OPE form of the B-meson LCD A was derived as [2J 



~t 2 A 2 



a s CF 



4tt 



2L 2 + 2L 



4A 
- it — 
3 



1 - 



a,C F { 2 16 35 
3 9 



a s C G 

47T 



-in 
4 2 



6 



1 - 



5^ 

a s CF 

47T 



o. s Cf 

47T 

t 2 A|(/i) 



22/ 



■4L-H 
4 



21/ 



4-7T 

l2~ 



5tt^ 
12 



2L 1 

3 12 



8tt 



(L-l) 



(65) 



with a s = a s (/i), as usual. This OPE form enables us to evaluate the J3-meson LCDA for interquark distances t 
with t < in a rigorous way in terms of three nonpcrturbative parameters in the HQET, i.e., a familiar HQET 
parameter, as the mass difference between the B-meson and &-quark, 



A = niB — TUb , 



(66) 



which is known to be associated with matrix elements of dimension-4 operators [24|, and the novel HQET parame- 
ters A|;(/x) and X^r(u), which are defined by matrix elements of the quark-antiquark-gluon three-body operators of 
dimension 5 as [13, un, ^3 



(0\qgE ■ a. l5 h v \B{v)) = ^(m)A|(m) 



(0\qgH ■ cr l5 h v \B{v)) = iF(n)X 2 H {n) 



(67) 



associated with the chromoelectric and chromomagnetic fields, respectively, in the rest frame (v = (1,0)). We 
note that the "universal" double logarithmic term, — (oi s Cf /47r)2L 2 , in ([6"5]) yields as a UV-finite term from the 
contribution of the diagram of Fig. [T] (a) in the one-loop matching calculation of the Wilson coefficients, while the 
UV-divergent part from the same diagram induced the cusp anomalous dimension in the evolution equation ([5)1 
through the renormalization constant in (j?]), as discussed in Sec. [TTJ Indeed, the derivative of the double logarithmic 
terms in ([65]) with respect to fi reproduces the term associated with the cusp anomalous dimension in fl8j). Taking also 
the derivative of the other terms in ([6"5"]) and combining the result with the scale dependence of the HQET parameters 
of ([66]) and ([67)) . i.e., dA/dfi = and that of X 2 ^ #(m) in terms of the the one- loop mixing matrix obtained in [35| . 
one can show that ([65)) satisfies the evolution equation ([8]), as demonstrated in (23[. Namely, the OPE form ([65]) 
for the LCDA is completely consistent with the RG structure of ([T]) embodied in ([5]). This important property can 
be demonstrated in an alternative way using our NLL-level solution (|35[) with (I43p -(I46 [) : expanding this solution in 
powers of a s (n), we obtain 



4>+(-ir, n) 



a s C F . 2 M a sC F , M 
1 In In — 



2tt 



Mo 

i 



2tt 



/'o 



^— - ln(r^ e 7E ) In — 

7T fl 



4>+{-ir, Mo) 



a s C F ln _M_ i dz(t> + (-~iTz,no) 
7T Mo Jo V 1 



(68) 



where the anomalous dimensions ©, ([TU]) are substituted, and the ellipses stand for the terms of order a 2 , and higher. 
Calculating the RHS of ([55]) with the substitutions, <j>+(— ir,Mo) — > 0+ PE ( — Mo): 4>+{—i TZ , Mo) ~> 0° PE (— irz, /io), 
we find that the result to 0(a s ) reproduces exactly </>+ PE ( — m) given by ([65]) . 

From the above discussions, our formula ([35]) with ([43 ]) -([46 ]) and with the input DA as <j>+(— irz, Mo) — ^ 
^+ PE (— *tz, /ito) embodies the model-independent form of B-meson LCDA |T]), incorporating the RG improvement 
associated with the NLL-level resummation, and satisfying that its expansion to the fixed order a s coincides with 
([6"5"]) . We emphasize that this form for the LCDA at the scale ijl is useful when r < 1/mo : recall that the OPE ([6"5"]) 
with the local operators and the corresponding Wilson coefficients calculated in the fixed-order perturbation theory is 
useful when the interquark separation is less than the typical distance scale of quantum fluctuation, i.e., r < 1/m [23J; 
but, the quasilocal structure in (|35p . manifested as its integral over < z < 1, reveals that this formula is applicable 
when the OPE form ^° PE ( — * T2 i Mo) substituted into the input DA is useful. Namely, the range of applicability of 
the OPE (|65[) obtained in perturbation theory is extended by performing the RG resummation as l|35[) to all orders 
in a s . In the following, we choose the input scale as /iq = 1 GeV, corresponding to the typical hadronic scale, and 
evaluate ([35]) with /i = /ihc ~ \/ m b^QCB- 

We now specify the values of the HQET parameters participating in our input DA, ([6"5]) with /i = /Lto- For those 
parameters, we adopt the same input values as used in our previous work [23j |: A of (|66p . which is defined by the 6-quark 
pole mass m^, is eliminated in favor of a short-distance parameter, Ada, free from IR renormalon ambiguities (36j 
and written as 



A = Ada(m) 



1 + 



7J3foJjj) 
16tt 



M" 



9CfQ: s (/z) 
"8^ 



(69) 




FIG. 2: The evolution of the B-meson LCDA based on the coordinate-space formula ()35[) : the solid curve shows the input 
DA, given by (a) 02 PE (— ir, /io) of (|65p and (b) its extension using (|71|l . The dashed, dot-dashed, and dotted curves show the 
results evolved to /i = 2.5 GeV, using the evolution operator of (|35|) with the NLL accuracy, using the evolution operator (|36[) 
corresponding to K = in (|17p . and using the LL-level evolution, respectively. 



to one-loop accuracy. Here, Ada(m) can be related to another short-distance mass parameter whose value is extracted 
from analysis of the spectra in inclusive decays B — >• X s ^y and B — >• A U Z i/, leading to Ada(mo) — 0.52 GeV. For the 
other two HQET parameters, we use the central values of 

A|(Mo) = 0.11 ±0.06 GeV 2 , A|(/i ) = 0.18 ± 0.07 GeV 2 , (70) 

which were obtained by QCD sum rules [HJ; at present, no other estimate exists for \ 2 E or \ 2 H . We now calculate 
(|65[) with (i = [io and the imaginary light-cone separation as t — > —it, and obtain model-independent description of 
the B-meson LCDA (j)+ FE (— it, (iq) for r < l/fJ*o = 1 GeV -1 , which is displayed by the solid line in Fig. [2] (a). This 
result can be substituted directly into the RHS of ([55]) as the input LCDA for the case with r < 1 GeV -1 , because 
zt < t < 1 GeV -1 = l//io in the integrand, reflecting the quasilocal nature as noted above. 

We now discuss the results of our evolution (J35J) to higher scale /i ~ ^/ 'nib Aqcd, shown in Fig. [2] (a): the dashed line 
denotes the full result of the LCDA </> + (— it, ^) at u, = 2.5 GeV, obtained by our NLL evolution (|55"]) using P5|) - (|ir]| . 
When we omit the effect of the DGLAP-type kernel K(z,a s ) in the evolution equation (|T7|). the resultant evolution 
is induced only by the factor (see (|5T)) ). yielding the dot-dashed curve. Omitting the other NLL terms in ([55)1 
further, as V(/i, /io) - ► (47r/ci! s (/z))/i( )(x) and £ — ► 0, we obtain the dotted curve that corresponds to the result of the 
LL-lcvel evolution. We see the considerable Sudakov suppression not only at the LL level but also at the NLL level; 
in particular, at the NLL level, the suppression arises in the moderate r regions, while the DA is enhanced for small 
r, reflecting the r dependence of the factor (|36|) . On the other hand, the DGLAP-type kernel contributes to shifting 
the distribution from small to moderate r, as a result of the integral over z in (|35| : such effect is characteristic of 
the evolution that is induced by the kernel associated with the plus-distribution of the type ([5]), and is similar to the 
corresponding effects arising in the usual DGLAP equation for the parton distribution functions of the nucleon. From 
the discussion above (|69[) . our full result, the dashed curve, is useful for providing model-independent behavior of the 
B-meson LCDA in small and moderate r regions presented in Fig. [2] (a), where the solid curve using the OPE form 
4>+ PE (— i T >Mo) is suitable for the input DA. 

The OPE form ([63]) . used for the input DA, breaks down in the large r region, where the contributions associated 
with the operators of any higher dimension become important because the contributions from the dimension-d oper- 
ators grow as ~ T d " 3 with increasing r. According to our previous work p3l |. we rely on a model function to describe 
the DA in the large r region dominated by the nonperturbativc effects, and, specifically, we use the following form of 
the input DA for the entire range of r, 

M-ir, Mo ) = 6(t c - r)0g PE (- t r, Mo ) + 6(r - r c ) ^ , (71) 

with t c ~ 1/no, such that we connect the small and moderate r behavior given by the rigorous OPE form, the first 
term, smoothly to a model for the large r behavior, the second term. This second term corresponds to the exponential 
form N(oj /wo)e~ u ^° in the momentum representation, and such form was suggested in an estimate of the i?-meson 
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A s i ( M ) [GeV" 1 ] 


H [GeV] 


Eq. H72J with Eqs. J35J, ([71| 


Lee-Neubert 


Braun et al. 


1.0 
1.5 
2.0 
2.5 


2.7 (0.6 + 2.1) 
2.4 (0.6 + 1.8) 
2.2 (0.5 + 1.7) 
2.1 (0.5 + 1.6) 


2.1 
1.9 
1.7 
1.6 


2.2 
2.0 
1.9 
1.8 



TABLE I: The results of the inverse moment A^ (fi) using (|35|) to the NLL accuracy for the input DA (|71[) . with the first and 
second numbers in the parentheses denoting the contributions from the first and the second terms in ((72} ■ The results obtained 
by Lee and Neubert [20l | and the estimates based on the calculation by Braun et al. [lf| are also shown for comparison. 

LCDA using QCD sum rules [l2| and was also adopted in [2(J as a nonperturbative component to model the i?-meson 
LCDA using the information from the OPE with the local operators of dimension d < 4 and the NLO corrections to the 
corresponding Wilson coefficients taken into account. (For the correspondence and difference between our OPE (p5|) 
and the OPE derived in [2(|, see the discussion in (Ij|.) Here, the two parameters N and ujq are determined by the 
continuity of (|71[) and its derivative, 4>+{— it, fio) and <90+(— it, /xo)/9t, at r = t c . The resulting values N ~ 0.86 and 
u>o ~ 0.31 GeV are found to be stable under the variation of r c for 0.6 GcV -1 < r c < 1 GcV -1 , and so is the behavior 
of the corresponding DA (|7Tj) (23|. In the following, we take r c = 1 GeV -1 , and now the solid curve in Fig. [5] (a) 
is continued to the r > 1 GeV -1 region with (|71l) . as presented by the solid curve in Fig. [2] (b). Using this result 
of (|7Tj) as the input DA in (j3"5j) , we obtain the other curves in Fig. [2] (b), which are evolved in the same way as the 
corresponding curves in Fig. [5] (a); in particular, the behaviors of those new curves in the region r < t c = 1 GeV -1 
completely coincide with those of the corresponding curves in Fig. [5] (a). Namely, the model-independent nature for 
t < t c , originating from the OPE, is preserved under the evolution. This remarkable feature of our results is a direct 
consequence of the quasilocal structure of the evolution (|35[) in the coordinate-space representation: the results in 
the region t < t c are not contaminated under the evolution by the assumed model behavior for larger distances, the 
contribution of the second term of (|7ip . On the other hand, the smooth continuation of this second term of (f7T]) to 
the first term at r = r c and the evolutions of the result lead to the similar interrelations between the curves at large 
t > t c in Fig. [21 (b) as those at moderate r, displayed also in Fig. [2] (a); e.g., we observe the considerable Sudakov 
suppression also at large r region. 

Integrating the dashed curve in Fig. [2] (b) over the entire range of r, we obtain the first inverse moment of the 
S-meson LCDA (see (07])-®), 

/•Tc _ POO 

X b 1 (l l ) = J dT( t>+(- iT ^) + J dT(j) + (-iT,n) , (72) 

at the scale jj, = 2.5 GeV, as the sum of the model-independent contribution, the first term, originating from the 
OPE and the second term depending on the assumed model behavior at large distances. Table Q] shows the results 
for fx = 2.5 GeV and some other values of li, with the first and second numbers in the parentheses denoting the 
contributions from the first and second terms in ([72")) . One can check that the exactly same values of A^ 1 (//) are 
obtained using the formula (|53[) . as the integrals of the input DA (|7ip with the corresponding weight function at 
the NLL accuracy. In Table [H we note that the result for ii = /j,q = 1.0 GeV coincides with that reported in our 
previous work [23. The evolution decreases A^ 1 (^) with increasing li, in particular, through the decrease of the 
model-dependent contribution, the second term of (I72|) . On the other hand, our results of A^ (li) for \i ~ m;,AQCD 
are larger than the results of [2(3] as well as of |l6|. where, for the former case in Table U we quote the results 
calculated in [2(j, and, for the latter case, we present estimates with the fixed-order formula ([FIT]) substituting A^ 1 (^t ) 
and <7i (/xo) (see (|49| ) that were obtained in [lf|. We recognize that the evolution from /i = fj,o to /^hc ~ \/ttmJ^qcd 
could give rise to the decrease of A^(/z) by 20-30%, with the larger value of A^ 1 (/x ) leading to the larger A^ 1 (/Lth c ); as 
emphasized in [23], our larger A7 3 1 (/j, ) than the corresponding values of other works [jillioj] originates from the novel 
contribution of Ag and \ 2 H in the OPE form (|65[) . which are associated with the dimension-5 operators representing 
the quark-antiquark-gluon three-body correlation (sec (J57J)). 

Our results of A7 J 1 (/x) presented in Table Q] as well as those calculated at even higher /z are plotted by the solid curve 
in the first panel in Fig. [3j and, similarly, the solid curves in the other two panels show the behaviors of the logarithmic 
moments defined in the coordinate space, R\{li) and i?2(/^) of (f54|) and (|55|) . with the NLL accuracy (j43 |) -(|46 |) using 
the input DA (fTTj). Here, the dot-dashed curves present the fixed-order calculations based on (|ol"j) - (|fJ3")) substituting 
Ag 1 (^o) I Ri((J>o), ^2(^0) and i?3(//o) calculated with (fTTj) . and these results are modified into the dashed curves when 
we omit the double logarithmic correction behaving as oc a s ln 2 (/z//io), compared to the corresponding tree (O(o^)) 
contribution, in each coefficient of A^ 1 (/x ), Ri{liq), i?2(Mo) and i?3(/io) in the RHS of ([pTj) - (|55|) . Furthermore, those 




H [GeV] |i [GeV] 

FIG. 3: The evolution of the first few logarithmic moments defined in the coordinate space, f|48f) . using (|71|l as the input 
B-meson LCD A at the scale /io = 1 GeV: the solid curves show the results of (|53|) -(|55 |) at the NLL accuracy, the dot-dashed 
curves show the fixed-order results (|61[) - (|63|1 . and the dotted curves show the results when a s — > 0. The dashed curves show 
the fixed-order results (f6lj) - (|63)l with the dominant double-logarithmic corrections omitted. 

results reduce to the dotted lines when a s — > 0. We find good accuracy of the fixed-order formulae (| 6 II) - (163)) . and thus 
the rapid convergence of the resummed perturbation theory in with when organized by \ of (|4ip . 

whose behavior as a function of fi is shown by the solid curve in Fig. |4] On the other hand, in ([6T]) -([63 j) . the double 
logarithmic effects play important roles to determine the scale dependence of Ag 1 (/z), iJaf/i), while the other 

0(a s ) contributions tend to cancel to a large extent. As the result, the perturbativc evolution from /i = /io to /% c 
can modify the values of those logarithmic moments considerably, by 20-30%. In Fig. [4j we also show the behavior 
of £ of (|46|) by the dashed curve; this demonstrates that the condition ([60]) is indeed satisfied, so that our formulae 
(f53")) - (|55p . as well as ([33)1 giving their basis, describe the well-defined evolutions for all relevant scales. 

In the present paper, we have discussed in detail the effects of the evolution on the S-mcson LCDA and on its 
integrals X^ 1 (fi), Ri(ii), i?2(A0 relevant to exclusive B decays, emphasizing modcl-indcpcndent aspects revealed by 
our coordinate-space approach, but did not intend to determine the precise values of those integrals at /i = /%<;. Such 
determination requires us to calculate the B- meson LCDA ((T|) at the initial scale reducing the corresponding 
theoretical uncertainty as much as possible: as found in our previous work [23[, the corresponding initial LCDA, 
calculated in the form of (fTTj) . is significantly influenced by the novel HQET parameters Xe and Ah arising in the 
OPE form (|65| . which arc associated with the dimension-5 quark-antiquark-gluon operators. Therefore, the rather 
large uncertainty in their existing estimate (|70[) based on QCD sum rules calls for more precise estimates of A# and 
Ah- Recently, higher-order corrections to the QCD sum rules for and Ah have been calculated, and these new 




H [GeV] 



FIG. 4: x °f l|4ip and £ of (|46[) as functions of /i with /io = 1 GeV, shown by the solid and dashed curves, respectively. 



contributions arc found to improve the estimate of A^ and A# [37[ . Wc also note that in the RHS of (|72j) , evaluated in 
TablcUl the second term is much larger than the first term. This suggests that As, as well as -Ri,2, is rather sensitive 
to the functional form that models the LCDA in the long-distance region; for example, a functional form motivated 
by the so-called Wandzura-Wilczek approximation 2M> 1221 provides an interesting possible alternative to the form 
appearing in the second term in (fTTj) (see, e.g., [38[ for other studies on the behaviors of the LCDA). Systematic 
investigations of these points, combined with the evolution effects obtained in this paper, could determine the values 
of A^ 1 (//), Ri(fj,), i?2(^) &t fi = /ihc as precisely as possible, and those results will be presented elsewhere. 



VI. CONCLUSIONS 



In this paper, we have studied the RG evolution of the B-meson LCDA, working in the coordinate-space represen- 
tation of the LCDA. The corresponding evolution equation and its solution demonstrated that our coordinate-space 
approach has remarkable advantages over the conventional approach in the momentum space. Indeed, only in the co- 
ordinate space, the relevant kernel in the evolution equation, associated with the cusp anomalous dimension as well as 
the DGLAP-typc anomalous dimension, is quasilocal, and this quasilocality is inherited by the corresponding analytic 
solution, leading to the simplest expression possible for calculating the evolution of the B-meson LCDA. Our explicit 
formula of the solution has the accuracy beyond the one-loop level in the RG-improved perturbation theory, taking 
into account the effect of the two-loop cusp anomalous dimension according to consistent order counting, such that 
the Sudakov-type double logarithmic effects as well as the DGLAP-type single-logarithmic corrections are resummed 
at the NLL accuracy. This result, in turn, allowed us to derive the master formula, by which the relevant integrals 
of the LCDA at the scale [ihc ~ \nnb^QCU, arising in the factorization formula for the exclusive B-meson decays, 
can be reexpressed in a model-independent way by the compact integrals of the LCDA at a typical hadronic scale 
fi ~ 1 GeV. 

We applied our evolution formula to the LCDA with the initial scale [1$, which is determined by the OPE having the 
perturbative (NLO) accuracy consistent with the NLL-level resummation, and the highest nonperturbative accuracy, 
at present, taking into account the local operators of dimension d < 5. The quasilocal structure of our evolution 
guarantees that the B-meson LCDA at a certain quark-antiquark distance is not contaminated under the change of 
the renormalization scale by the configurations of the quark and antiquark for the larger distances, so that the LCDA 
at high scales, obtained through our evolution from the OPE-based, initial LCDA that is accurate for interquark 
distances less than ~ l/fi ~ 1 GeV -1 , exhibits the model-independent behaviors for distances < 1 GeV -1 . Our 
explicit numerical calculation indicated the considerable effects of the evolution, from the scale fig to fihe, for the 
LCDA and its relevant integrals. In particular, we found that the dominant roles are played by the double logarithmic 
corrections, although particular attention was not paid to them in previous works. On the other hand, we observed the 
rapid convergence of the corresponding resummed perturbation series organized by the proper logarithmic expansion. 

Using the information available for the nonperturbative effects associated with the OPE-based, initial LCDA, 
our evolution gave an estimate for the relevant integrals of the LCDA at the scale /ihc, e.g., Ag 1 (/i) ~ 2.1 GeV -1 at 
fj, = 2.5 GeV. This is larger than the estimates by other works, inheriting the larger value A^ 1 (fi = 1 GeV) ~ 2.7 GeV -1 
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in our case, which is induced by matrix elements of the dimension-5 quark-antiquark-gluon operators in the OPE for 
the initial LCDA. Combined with an update of the information on the nonperturbative effects in the initial LCDA, the 
results derived in this paper are immediately applicable for calculating the refined values of those integrals relevant 
to exclusive B decays. 



Appendix A: The evolution in the momentum representation 

The evolution of the i?-meson LCDA </>+(i, (i) in the representation with the light-cone separation t is provided by 
our solution (|35[) with the analytic continuation r — > i(t — iO) performed. We calculate the Fourier transformation of 
this result, in order to derive the evolution for the LCDA </>+(w, /i) in the momentum representation (see (U}): 

If 00 _ 6 V(ai^o) + (1-27e)? r°° 

4> + {u,n) = — die w V+(t,M) = - — -z / du'Ifaw'W+iu',!*), (Al) 

27r J-oo r(£)/4 Jo 

with the integration kernel, 

I(u,u')= dz(-^—\ / — e^-^^t-iO)] - * , (A2) 



J-zJ J-oo 2tt 
where the integration over t can be performed straightforwardly, yielding 

1 f 1 / z 0(lu-lu'z) 



1 



m J 



dzz 1 -^! - z)*- l (w - u'zf- 1 . (A3) 



Here, we have introduced the notation, uj < = min(a;,a/). Changing the integration variable to u = {u)'/u)<)z and 
using louj' = uj > uj < with oj> = max(w,u'), we can rewrite (|A3|) as 



IM = £^ fduu^Hl-uf-'il-^u)'' 1 , (A4) 

r(0 Jo \ u> j 

and we note that this can be expressed by the hypergeomctric function, 

F , a , F( 7 ) ^ r(a + n)r(P + n) z n 

2 Fi a,/3; 7 ;z ) = — — — — - > — — ■ — A5 

r (a) r (/3) ^ r( 7 + n) n\ 



I" duu fj - 1 (l- U y- fj - 1 (l-zu)- a , (A6) 



r(/?)r( 7 - p) J 



where the first line shows the usual definition by the series expansion, and the second line gives the integral represen- 
tation to be compared with (|A4|) . Substituting the result into (|Alj) . we obtain 



<M«, M) = e^^-^^-f ^ ^ + (u/, Mo) ^ a *k f 1 - e, 2 - f ; 2; ^ ) , (A7) 

which gives a well-defined formula when the condition (|3"Tj) is satisfied. The evolution of the B-meson LCDA in this 
form was first derived in [H|,[2(| by solving the evolution equation given in the momentum space, see (TT3")) , (fT4"|); note 
that V(/i,/io) +£ and £ in the present paper correspond, respectively, to V(fj,,/io) and g in [18|, |20(. We mention that 
it would not be straightforward to derive the relations (|56l) -([59l) based on (|A7[) . because of the structure involving the 
complicated integration of the hypcrgcometric function (see [34|V 
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